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We study bit-parallel-wavelength (BPW) pulse transmission in multi-channel single-mode optical 
fiber links for high-performance computer networks. We develop a theory of the pulse shepherding 
effect earlier discovered in numerical simulations, and also describe families of the BPW solitons 
and bifurcation cascades in a system of N coupled nonlinear Schrodinger equations. 



A growing demand for high-speed computer communi- 
cations requires an effective and inexpensive parallel com- 
puter interconnects that eliminate bottlenecks caused 
by parallel-to-serial conversion. Recently, bit-parallcl- 
wavelength (BPW) single-fiber optical links were pro- 
posed as possible high-speed computer interconnects . 
Functionally, the BPW link is just a single-medium par- 
allel fiber ribbon cable. In a conventional ribbon cable, N 
parallel bits coming from a computer bus are transmitted 
by N pulses travelling in N separate fibers. In a BPW 
scheme, all N bits are wavelength multiplexed and trans- 
mitted by time-aligned pulses, ideally - solitons, through 
a single optical fiber. 

For any bit-parallel transmission the crucial problem 
is mantaining the alignment of pulses corresponding to 
parallel bits of the same word. Unlike the fiber ribbon ca- 
ble, a single-fiber BPW link presents a unique possibility 
of dynamical control of the pulse alignment by employ- 
ing the so-called pulse shepherding effect when a 
strong "shepherd" pulse enables manipulation and con- 
trol of co-propagating weaker pulses. Experimentally, the 
reduction of the misalignment due to the group-velocity 
mismatch of two pulses in the presence of the shepherding 
pulse has been observed in a dispersion-shifted Corning 
fiber j§. 

In this Letter, we develop a rigorous theory of the shep- 
herding effect, and show that it is caused by the nonlin- 
ear cross-phase modulation (XPM) of pulses transmitted 
through the same fiber. The co-propagating pulses can be 
treated as fundamental modes of different colour trapped 
and guided by an effective waveguide induced by the the 
strong shepherd pulse. The resulting multi- component 
soliton pulse propagates in the fiber preserving the time 
alignment of its constituents and thus enables multi- 
channel bit-parallel transmission. For the first time to 
our knowledge, we analyze multi-component solitons and 
describe a mechanism of the soliton multiplication via the 
bifurcation cascades in a model of N nonintegrable cou- 
pled nonlinear Schrodinger (NLS) equations. 

To describe the simultaneous transmission of N pulses 
of different wavelengths in a single-mode optical fiber, we 
follow the standard derivation j^] and obtain a system 
of N coupled NLS equations in the coordinate system 
moving with the group velocity v g o of the central pulse 



{o<j<N-iy, 

where S m j = 2 I]m#j(7m/7j)|A™| 2 . For a pulse j, 
Aj(z,t) is the slowly varying envelope measured in the 
units of \fPo, where Pq is the incident power carried 
by the central pulse, aj = — /?2j / 1 /?20 1 is the normal- 
ized group velocity dispersion, 8j = (v g o — v gj)/v g ov g j 
is the relative group velocity mismatch, and jj — uJj/uiQ 
characterises the nonlinearity strength (cvo = 70 = 1). 
The time, t = (T — Z/v g o)/T , and propagation dis- 
tance, z = Z/Lq, are measured in units of character- 
istic pulse width, To « 10 ps, and dispersion length 
L o = Tq/\P 2 o\ ~ 50 km [HQ. For the operating wave- 
lengths spaced 4-^5 nm apart (to avoid the four-wave- 
mixing effect), within the band 1530-^1560 nm (see Refs. 
HH), the coefficients aj and jj are different but close to 
1 . For the realistic group- velocity difference of less than 
5 ps/km Q], the mismatch parameter 8j < 1. 

Below we show that the system (|l|) admits station- 
ary solutions in the form of multi-component BPW soli- 
tons which represent a shepherd pulse time-aligned with 
a number of lower-amplitude pulses. We then discuss the 
effect of the group-velocity mismatch on the time align- 
ment of the constituents of the multi-component pulse. 

To find the stationary solutions of Eqs. ([j]), we use the 
transformation: Aj(t, z) = Uj(t)exp(iSjajt + iXjZ), with 
the amplitudes Uj obeying the following equations: 

1 d 2 uo j 2 21 1 

2~aW + I U ° + ^ lnUn j U ° = 2 U °' 

, 2 / N-l \ 

where the amplitudes, time, and A n are measured in units 
of \/2Ao, (2Ao) -1 / 2 , and 2Ao, respectively. 

System (^|) has exact analytical solutions for N cou- 
pled components. Indeed, looking for solutions in the 
form uo(t) = C/osechi, u n (t) = t/ ra sech t, we obtain the 
constraints A„ = a n /2, and a system of N coupled alge- 
braic equations, 
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Ul + 2 ]T InUl = 1, 7n [/ 2 + 2 £ 7m C/2 = a n . 

n—l m^n 

As long as all modal parameters, a n , 7n , are close to 1, 
this solution describes a composite pulse with N nearly 
equal constituents. In the degenerate case, «„ = 7n = 1, 
the amplitudes are §: f7 = U n = [1 + 2{N - 1)]~ 1/2 . 

Approximate analytical solutions of different types can 
be obtained in the linear limit 0, when the central fre- 
quency pulse, uq, is much larger than other pulses, and 
the XPM interaction between the latter can be neglected. 
Lineariziation of Eqs. (Q) for \u n \ -C |uo| yields an exactly 
solvable NLS equation, for uo, and N—l decoupled linear 
equations, for u n . Each of the latter possesses a localized 
solution provided A„ = A„ = (a„/8)[l — \/l + 16a„J 2 . 
Near this point, the central soliton pulse (shepherd pulse) 
can be thought of as inducing an effective waveguide that 
supports a fundamental mode u n , with the corresponding 
cutoff, A„. Since a n and 7n are close to 1, the soliton- 
induced waveguide always supports no more than two 
modes of the same wavelength, with largely separated 
eigenvalues. As a result, the effective waveguide induced 
by the shepherd pulse stays predominantly single-moded 
for all operating wavelengths. 

Let us describe the mechanism of pulse shepherd- 
ing in more details. First, we consider the simplest 
case N = 2. If the two pulses, (0) and (1), do not 
interact, then the uncoupled Eqs. (Q) possess only 
single-component soliton solutions Uo(t) — sechi and 
ui = (2Ai 7 i/ai) 1 / 2 sechy / 2Ai/oii, with the correspond- 
ing normalized powers Pq = J u^dt = 2 and Pi = 

2(2Ai 7 i) 1 / 2 a^ 3/2 . These solutions can be characterised 
by the curves on the diagram P(X) [see curves (0) and 
(1) in Fig. 1]. When the two copropagating pulses in- 
teract, a new branch of the two-mode solitons (0 + 1) 
appears (branch A-B in Fig. 1). It is characterised by 
the total power P(Ai) = Pq + Pi, and it joins the two 
branches -Po(Ai) and Pi(Ai) at the bifurcation points Oi 
and O2, respectively. Near the point Oi, the solution 
consists of a large pulse of the central wavelength that 
guides a small component u\ (see Fig. 1, point A). The 
point Oi therefore coincides with the cutoff Ai = Ai for 
the fundamental mode ui guided by the shepherd pulse 
uq. Shapes and amplitudes of the soliton components 
evolve with changing Ai (see the point B in Fig. 1). The 
component Mo disappears at the bifurcation point Oi- 

Next, we consider a shepherd pulse guiding three 
pulses, n = 1,2,3, with the modal parameters: 00,3 = 
70,3 = 1, ai = 7 i = 0.65, and ai = 7 2 = 0.81. Soli- 
tary waves of this four-mode BPW system can be found 
as localized solutions of Eqs. (^), numerically. Figure 2 
presents the lowest-order families of such localized solu- 
tions on the line Ai = A2 = A3 = A in the parameter space 
{Ai,A2,A3}. If the pulses (1), (2), and (3), were interact- 
ing only with the central pulse but not with each other, 



then the bifurcation pattern for each of this pulses would 
be similar to that shown in Fig 1. Thin dotted, dashed, 
and dash-dotted curves in Fig. 2 correspond to the soli- 
tons of three independent pulses (1), (2), and (3), shown 
with branches of corresponding two-mode solitons of the 
BPW system with pairwise interactions, (0+1), (0 + 2), 
and (0 + 3), respectively (cf. Fig. 1). In fact, all four 
pulses interact with each other, and therefore each new 
constituent added to a multi-component pulse "feels" the 
effective potential formed not only by the shepherd pulse 
but also by all the weaker pulses that are already trapped. 
In addition, mutual trapping of the pulses (1), (2) and (3) 
without the shepherd pulse is possible. As a result, new 
families of the two-mode (1 + 2) and, branching off from 
it, the three-mode (0+1 + 2) solutions appear (marked 
curves in Fig. 2). The three- mode solutions bifurcate at 
the point O3 and give birth to the four-mode (0+1+2+3) 
solitons (branch O3 — O4). An example of such four- wave 
composite solitons is shown in Fig. 2 (inset). This solu- 
tion corresponds to the typical shepherding regime of the 
BPW transmission for N = 4, when the central pulse uo 
traps and guides three smaller fundamental-mode pulses 
on different wavelengths. 

On the bifurcation diagram (Fig. 2), starting from 
the central pulse branch, the solution family undergoes a 
cascade of bifurcations: 0\ — > Oi O3 —>■ O4. On each 
segment of the corresponding solution branches, differ- 
ent multi-component pulses are found: (0) — > (0 + 1) — > 
(0 + 1 + 2) -> (0 + 1 + 2 + 3) -► (1 + 2 + 3). The values 
of the modal parameters in Fig. 2 are chosen to provide 
a clear bifurcation picture, although they correspond to 
the wavelength spacing that is much larger than the one 
used in the experiments for which 7n / 7 n+i ~ 0.997. 
If the modal parameters are tuned closer to each other, 
the first two links of the bifurcation cascade tend to dis- 
appear. Ultimately, for equal parameters, the bifurcation 
points O2 and O3 merge at the point Oi, and the four- 
mode soliton family (thick line in Fig 2) branches off di- 
rectly from the central pulse branch (0). Then, near the 
point Oi, the four- mode pulse can be described by the 
linear theory. The qualitative picture of the bifurcation 
cascade for N — 4 preserves for other values of N. 

The BPW solitons supported by shepherd pulses are 
linearly stable ||. However, the effect of the relative 
walk-off due to the group velocity mismatch of the soli- 
ton constituents endangers the integrity of a composite 
soliton and thus the pulse alignment in the BPW links. 
It is known that, in the case of two solitons of comparable 
amplitude, nonlinearity can provide an effective trapping 
mechanism to keep the pulses together ||. In the shep- 
herding regime, a multi-component pulse creates an ef- 
fective attractive potential, and the j-th pulse is trapped 
if its group velocity is less than the escape velocity of this 
potential. The threshold value of the walk-off parameter 
can be estimated as f|: (5 2 < (4a>j/Pj) J2 m & 7m«m(°)> 
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where u, n (0) is the peak intensity of the component m. 
For instance, for the component j = 1 of the four- 
component BPW soliton presented in Fig 2 (point A), 
the estimated threshold <5i < 1.7 agrees with the numer- 
ically calculated value Si < 2.2. 

In reality, all the components of a BPW soliton would 
have nonzero walk-off. The corresponding numerical sim- 
ulations are presented in Fig. 3 for TV = 4. Initially, we 
launch four pulses as an exact four-mode BPW soliton 
A (see Fig. 2) of Eqs. (||). When this soliton evolves 
along the fiber length, z, in the presence of small to 
moderate relative walk-off (Sj ^ for j + 1 0), its com- 
ponents remain strongly localized and mutually trapped 
[Figs. 3(a,b)], whereas it loses more energy into radia- 
tion for much larger values of the relative walk-off [Figs. 
3(c,d)]. The former situation is more likely to be real- 
ized experimentally as the relative group-velocity mis- 
match for pulses of different wavelength is different [Q. 
In this case, the conclusive estimate for the threshold 
values of Sj can only be given if the shepherd pulse is 
much stronger than the guided pulses, which are approx- 
imately treated as non-interacting fundamental modes of 
the effective waveguide induced by the shepherd pulse. 

In conclusion, we have developed a theory of the shep- 
herding effect in BPW fiber links and established that the 
pulse shepherding can enable the time alignment of the 
copropagating pulses despite the relative group velocity 
mismatch. 
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of the Performance and Planning Fund. 




FIG. 1. Bifurcation diagram of the two-channel BPW 
model for qo,i = 7o,i = 1.0: horizontal line - branch of the 
central pulse (0); dashed - branches of the (1) and (0+1) 
solitons (A — B). Inset: examples of the (0+1) solitons cor- 
responding to the points A, B, shown for the uo (solid) and 
ui (dashed) components. 
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FIG. 2. Bifurcation diagram of the four-channel BPW 
model: horizontal line - branch of the central pulse (0); dot- 
ted, dashed and dash-dotted - branches of pulses (1), (2), and 
(3) with the (0+1), (0+2), and (0+3) solitons, respectively; 
solid thin - branches of the (1+2), (0+1+2), and (1+2+3) 
solitons; thick - branch of the (0+1+2+3) solitons. Inset: ex- 
ample of the (0+1+2+3) soliton corresponding to the point 
A; solid, dotted, dashed and dash-dotted lines - amplitudes 
of the uo, ui, v,2, and 113 components, respectively. 
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FIG. 3. Influence of the mode walk-off on the shepherding 
effect for the input parameters as for soliton A of Fig. 2. 
Shown are: components of the output pulse at 2 = 50 for (a) 
<5i = 0.45, 82 = -0.35, 5 3 = 0.25, and (c) 81 = S 2 = S 3 = 0.9; 
(b,d) - evolution of the position of the maxima for all pulse 
constituents. Solid line - maximum of the shepherd pulse. 
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